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Abstract 

By taking into account the contamination of foreground radiations, we employ the Fisher matrix 
to forecast the future sensitivity on the tilt of power spectrum of primordial tensor perturbations 
for several ground-based (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), balloon- 
borne (EBEX, Spider) and satellite (CMBPol, COrE, LiteBIRD) experiments of B-mode po¬ 
larizations. For the fiducial model nt = 0, our results show that the satellite experiments give 
good sensitivity on the tensor tilt nt to the level ant ^0-1 ^ ^ 2 x 10“^, while the ground- 

based and balloon-borne experiments give worse sensitivity. By considering the BICEP2/Keck 
Array and Planck (BKP) constraint on the tensor-to-scalar ratio r, we see that it is impossible 
for these experiments to test the consistency relation nt = —r/8 in the canonical single-field 
slow-roll inflation models. 
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1 Introduction 


The inflation paradigm E El El HI E] predicts the primordial gravitational waves, i.e. 
tensor perturbations. The primordial gravitational waves can contribute to the total 
intensity and polarizations of cosmic microwave background (CMB) anisotropy [SI El El 
El uni EH E2]- From the recent Planck TT,TE,EE+lowP data release [13], the upper 
limit on tensor-to-scalar ratio is given by ^ 0.002 < 0.11 at the 95% confidence level (C.L.) 
by fitting the ACDM+r model. The B-mode polarizations contributed by primordial 
gravitational waves may be detectable at the range ^ < 150. Actually, BICEP2 H has 
pushed the sensitivity of B-mode polarizations to be comparable with that of temperature 
in searching for the primordial gravitational waves. However, the polarized dust emissions 
make us difficult to distinguish whether the detected B-mode power comes from the 
primordial gravitational waves [IH EH EH EH]. Recently, Planck [TS] released the full-sky 
data of polarized dust emissions. Based on this, a joint analysis of the B-mode data from 
BICEP2/Keck Array and Planck (BKP) [20] yielded an upper bound ro .05 < 0.12 at the 
95% C.L., which is compatible with the upper limit from Planck data without the B-mode 
polarizations. 

The tilt of power spectrum of primordial tensor perturbations is used to measure the 
feature of the primordial gravitational waves. In the power spectrum of primordial tensor 
perturbations, the tensor tilt rit is defined by 

P,(k) = rA. (J-) ' , ( 1 ) 

where Pt{k) is the amplitude of power spectrum of primordial tensor perturbations at the 
scale k, r denotes the tensor-to-scalar ratio at a given pivot scale kp and Ag the amplitude 
of primordial scalar perturbations which is set as a constant in this paper. In the inflation 
model, the tensor tilt is generally predicted as nt = — 2e [211E2]. The inflation requires 
a/a = H‘^{1 — e) where e = H/H'^, and thus —2<nt< 0. In the canonical singe-field 
slow-roll inflation models, the tensor tilt is determined by the tensor-to-scalar ratio via 
the consistency relation, i.e. nt = —r/8 [21] • By considering the upper bounds on r, we 
expect the power spectrum of primordial gravitational waves to be nearly scale-invariant, 
i.e. nt ~ 0. The current constraint on nt was given by n^^o.oi = ~0.761 q;52 [23] at the 68% 
C.L. by combining only the BKP B-mode data and the upper limit on the intensity of 
stochastic gravitational wave background from Laser Interferometer Gravitational-Waves 
Observatory (LIGO) [23] • The scale-invariant spectrum is well compatible with this 
constraint. 

Even though at present there are no evidence for the primordial gravitational waves, 
several future polarization experiments might reach the sensitivity to detect the primordial 
gravitational waves in the coming years. As a recent analysis. Ref. I2SI forecasted that the 
primordial gravitational waves with theoretically motivated r ~ 2 x 10“^ can be achievable 
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by certain future experiments if the noise is reduced to ~ 1/ii^-arcmin and the lensing 
B-modes reduced to 10%. Their forecasts are not changed signihcantly with respect to 
previous estimates [26]. In this paper, we study the sensitivity on the tensor tilt rii for 
several future ground-based (AdvACT, CLASS, Keck/BICEP3, Simons Array, SPT-3G), 
balloon-borne (EBEX, Spider) and satellite (CMBPol, COrE, LiteBIRD) experiments. 
We are just interested in studying the nearly scale-invariant case, i.e. n* ~ 0 which 
corresponds to a class of the simplest inflation models. Similar to Ref. [25] , the approach 
of Fisher matrix is also used in our analysis. The paper is arranged as follows. The 
B-mode polarizations, foregrounds and noise sources are described in section]^ In section 
the method used in this paper is revealed. In section we show our forecasts for the 
future experiments. Conclusions and discussion are given in section 


2 Signal, Foregrounds and Noise 

In general, the CMB linear polarizations can be expressed in terms of the spin-weighted 
spherical harmonics ± 2 Wm, namely 


Q±iU = ±2Y£m , 

im 

or equivalently, dehned by the E- and B-modes as 

E = 'y ^ , 

im 

B = y ^ , 

where the coefficients are given by 


im 


^im 






1 




( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 


In this paper, we are just focused on studying the B-modes which may include the signal 
of primordial gravitational waves. 

In the linear perturbation theory, the primordial B-modes are Gaussian with zero mean, 
and their variance is given by 




4*m') = 


(7) 


where 6 comes from statistical isotropy. As conventions, the angular correlation coeffi¬ 
cients between B-modes are dehned as 


Cf = 


+ 7) ryBB 


2ti 


-a 




( 8 ) 
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where we have dropped the superscript BB for simplicity. In our study, we generate the 
power spectrum of primordial B-modes by running the CAMB package |27l|2H], and we 
set all the cosmological parameters except r and rit to the best-£t values of Planck 2015 
results [21]. However, we should find out a reasonable pivot scale kp such that there is 
least degeneracy between r and nt in Eq. ([^. 

The foregrounds contaminate the CMB B-mode signal and should be taken into account 
in the forecast. Actually, we can separate each component of foregrounds by noting that 
they have very different frequency-dependence. The Galactic synchrotron emission (S) 
and thermal dust emission (D) are considered in this paper. Their power spectra are given 
by 

Stu = h (9) 

Di, = (W°fc°^(W°fAD[^£Y , (10) 

where various parameters can be found in Tab. 0 and and are dehned by 


Parameters 

Synchrotron 

Thermal Dust 

^72% 

2.1 X 10“® 

0.169 

^53% 

2.1 X 10“® 

0.065 

^24% 

2.1 X 10-® 

0.019 

^11% 

4.2 X 10“® 

0.013 

^1% 

4.2 X 10"® 

0.006 

p[GHz] 

65 

353 

£ 

80 

80 

a 

-2.6 

-2.42 

/3 

-2.9 

1.59 

T[K] 

— 

19.6 


Table 1: A list of foreground parameters [161 [19], |25l EQ] . Here denotes the cleanest 
effective area fsky in the sky and its unit is 


= 

WCmb 

(c 




) ’ 

= 

WCmb 

I'D 

(C 

\ ^huo/ksT _ ^ 

WCMB 

\jyD, 

] ^hv/ksT _ 

CMB _ 



hv 

X = 7 - . 


(e® —1)^’ ksTcMB 


( 11 ) 

( 12 ) 

(13) 


Here we have rescaled the temperature of Galactic synchrotron emission and thermal 
dust emission with respect to the GMB temperature. There might be certain correlation 
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between synchrotron emission and dust emission. To account for this, we assume their 
correlation taking the form g in the power spectra. 

In this paper, we do not consider the systematics which strongly depend on the ex¬ 
perimental setups. However, we consider the instrumental white noise which is Gaussian. 
The power spectrum of the white noise can be expressed as [STj 


Af 


£(£+1) 

271 


6P 




(14) 


where 6P denotes the sensitivity for the Stokes parameters Q and U, and ab = 0A256fwhm 
denotes the beam-size variance. Various instrumental parameters here can be found in 
Tab. I^in Appendix. 

The gravitational lensing can also limit our ability to detect primordial B-mode polar¬ 
izations. However, the lensing B-modes and primordial B-modes have the same frequency- 
dependence, and we cannot separate them as mentioned above. Fortunately, one can re¬ 
construct the lensing potential by considering the CMB data of temperature and E-mode 
polarizations at small angular scales, and then remove the lensing B-modes away at large 
angular scales [42l 1431 HU |15]. In this paper, we assume the power of lensing B-modes 
reduced to 10% of its original value for the CMBPol and COrE experiments. For others, it 
is marginal to employ delensing [2S]- The residual power of lensing B-modes can 

be incorporated into the power spectrum of an effective noise, i.e. Me —> Me + 

m- 


3 Likelihood and Fisher matrix 


As mentioned above, we will deal with three components of CMB B-modes which 
originate from primordial gravitational waves. Galactic dust and synchrotron emissions, 
respectively. In the “Component Separation” (CS) method, the average of log-likelihood 
can be given by [25] 


(log Cbb) 


~2 ^)fsky 

e 


^log det 


fWCeW^ +Me\ 
\WCeW^ + Me) 


(WCeW^ + Me \\ 
\WCeWT + Me )) 


1 


( 16 ) 

where the bar denotes all the parameters hxed to their “true” values, and we have used 
the normalization such that (Yog Cbb) = 0 for Ce = Ce and W = W. Here W denotes the 
frequency-dependence of each component of CMB B-modes, which is a 3 x iV matrix with 
a row (1, hF,^). N denotes the number of frequency channels for each experiment. 

Ce is the covariance matrix of the amplitudes of three components, fsky stands for the 
effective area of the sky in Tab. where each experiment observes. 


In our consideration, the average of log-likelihood is a function of several parameters 
p which are {r,nt, AdMsMd, Ps, g)- We assume Gaussian priors for Ad, As, Pd and Ps 
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with the variance of 50%, 50%, 15% and 10%. Based on Eq. (15), the Fisher matrix is 
dehned by 

_ d^logCBB) I . . 

where p denote “true” values of the parameters p. In our hducial model, we set g = 0.5 
and fit = 0. The minimum error on the parameter pi is given by the Cramer-Rao bound, 

4 > ■ ( 17 ) 


In this paper, we will forecast the future sensitivity on both r and rit simultaneously, 
since primordial gravitational waves are not detected until now. There could be certain 
correlation between r and n*, or equivalently, the (r, rit) conhdence ellipse has a tilt. Thus 
we should hnd a pivot scale kp to make their correlation to be least. In other words, we 
should hnd the pivot scale such that = 0, which minimizes the constraint on r. 


4 Analysis and Results 

For the future ground-based experiments, we consider the CMB multipoles of the range 

[30.150] for Keck/BICEP3, Simons Array and SPT-3G, while [2,150] for AdvACT and 
CLASS. The reason is that AdvACT and CLASS cover a larger fraction of the sky. We 
forecast on the future sensitivity on r and rit by using the Fisher matrix. The la errors on 
r and rit can be found in Tab. We also list the pivot scale /Cp[Mpc“^] at which there is 
no tilt for the (r, rit) conhdence ellipse. We found that the ground-based experiments can 
probe the tensor-to-scalar ratio to the level r ~ 0.01, which is consistent with previous 
estimates. The la error on the tensor tilt is to the level a^ ~ 0.1 when 2 x 10“^ < r < 0.1 
for AdvACT and CLASS, while am ^ 1 for other three experiments. 

For the future balloon-borne experiments, we consider the CMB multipoles of the range 

[30.150] . The la errors on r and rit and the pivot scale kp can be found in Tab. We 
hnd that EBEX and Spider can marginally probe the tensor-to-scalar ratio to the level 
r ~ 0.02, which is consistent with previous estimates. The la error on the tensor tilt is 
am > 1 when 2 x 10“^ < r < 0.1. 

For the future satellite experiments, we consider the CMB multipoles of the range 
[2, 300] for CMBPol and COrE while [2,150] for LiteBIRD. The Icr errors on r and rit and 
the pivot scale kp can be found in Tab. The delensing has been taken into account for 
CMBPol and COrE. We hnd that the satellite experiments can well probe the tensor-to- 
scalar ratio even to the level r ~ 2 x 10“^, which is also consistent with previous estimates. 
The lex error on the tensor tilt rit runs from ~ 0.1 to ~ 0.01 when r runs from 2 x 10“^ 
to 0.1. Thus the future satellite experiments provide the highest sensitivity on the tensor 
tilt rit by contrast to the future ground-based and balloon-borne experiments. 
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r 


f^nt 

fcp[MpC 3] 


0.1 

5.3 X 10"^ 

9.7 X 10-2 

8.4 X 10-3 


0.05 

4.4 X 10-3 

1.2 X 10-3 

7.5 X 10-3 

AdvACT 

0.02 

3.6 X 10-3 

1.6 X 10-3 

5.2 X 10-3 

2.3 X 10-3 

0.01 

3.1 X 10-3 

3.4 X 10-3 


0.005 

2.6 X 10-3 

3.6 X 10-3 

2.1 X 10-3 


0.002 

2.1 X 10-3 

7.5 X 10-3 

1.3 X 10-3 


0.1 

6.7 X 10-3 

9.0 X 10-2 

6.3 X 10-3 


0.05 

5.5 X 10-3 

1.1 X 10-3 

4.7 X 10-3 

CLASS 

0.02 

4.1 X 10-3 

1.6 X 10-3 

2.4 X 10-3 

0.01 

3.1 X 10-3 

2.5 X 10-3 

1.3 X 10-3 


0.005 

2.2 X 10-3 

4.2 X 10-3 

8.6 X 10-3 


0.002 

1.5 X 10-3 

9.3 X 10-3 

5.7 X 10-3 


0.1 

2.3 X 10-^ 

1.2 

9.2 X 10-3 


0.05 

1.7 X 10-2 

1.8 

8.9 X 10-3 

Keck/BICEP3 

0.02 

1.4 X 10-2 

3.6 

8.7 X 10-3 

0.01 

1.3 X 10-2 

6.6 

8.7 X 10-3 


0.005 

1.2 X 10-2 

12.5 

8.6 X 10-3 


0.002 

1.2 X 10-2 

30.4 

8.6 X 10-3 


0.1 

1.7 X 10-2 

7.0 X 10-3 

9.3 X 10-3 


0.05 

1.5 X 10-2 

1.3 

9.3 X 10-3 

Simons Array 

0.02 

1.4 X 10-2 

3.0 

9.2 X 10-3 

0.01 

1.4 X 10-2 

5.9 

9.2 X 10-3 


0.005 

1.4 X 10-2 

11.7 

9.2 X 10-3 


0.002 

1.4 X 10-2 

29.1 

9.2 X 10-3 


0.1 

8.8 X 10-3 

4.3 X 10-3 

9.4 X 10-3 


0.05 

6.6 X 10-3 

6.4 X 10-3 

9.2 X 10-3 

SPT-3G 

0.02 

5.3 X 10-3 

1.3 

9.1 X 10-3 

0.01 

4.8 X 10-3 

2.3 

9.1 X 10-3 


0.005 

4.5 X 10-3 

4.4 

9.0 X 10-3 


0.002 

4.4 X 10-3 

10.6 

9.0 X 10-3 


Table 2: la errors on r and n* for future ground-based experiments. The pivot scale 
A;p[Mpc“^] is also listed, at which there is no correlation between r and n* and then the 
(r,nt) conhdence ellipse has no tilt. 
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r 


^nt 

fcp[Mpc 


0.1 

2.3 X 10 “^ 

1.2 

9.1 X 10 -^ 


0.05 

1.7 X 10-2 

1.9 

8.8 X 10-3 

EBEX 

0.02 

1.3 X 10-2 

3.7 

8.6 X 10-3 

0.01 

1.2 X 10-2 

6.8 

8.4 X 10-3 


0.005 

1.1 X 10-2 

12.9 

8.4 X 10-3 


0.002 

1.0 X 10-2 

31.4 

8.4 X 10-3 


0.1 

1.8 X 10-2 

1.0 

8.6 X 10-3 


0.05 

1.5 X 10-2 

1.8 

8.4 X 10-3 

Spider 

0.02 

1.4 X 10-2 

4.1 

8.3 X 10-3 

0.01 

1.3 X 10-2 

7.9 

8.3 X 10-3 


0.005 

1.3 X 10-2 

15.4 

8.3 X 10-3 


0.002 

1.3 X 10-2 

38.2 

8.3 X 10-3 


Table 3: la errors on r and n* for future balloon-borne experiments. The pivot scale 
A;p[Mpc“^] is also listed correspondingly. 



r 

'p 

<^nt 

fcp[MpC 3 ] 


0.1 

1.2 X 10-3 

3.0 X 10-2 

1.1 X 10-2 


0.05 

7.1 X 10 -^ 

3.8 X 10-2 

1.0 X 10-2 

CMBPol 

0.02 

3.8 X 10 -^ 

5.1 X 10-2 

9.0 X 10-3 

0.01 

2.5 X 10 -^ 

6.2 X 10-2 

8.4 X 10-3 


0.005 

1.8 X 10 -^ 

7.5 X 10-2 

7.9 X 10-3 


0.002 

1.4 X 10 -^ 

9.7 X 10-2 

7.2 X 10-3 


0.1 

1.3 X 10-3 

3.4 X 10-2 

1.1 X 10-2 


0.05 

8.2 X 10 -^ 

4.4 X 10-2 

9.7 X 10-3 

COrE 

0.02 

4.6 X 10 -^ 

5.8 X 10-2 

8.6 X 10-3 

0.01 

3.3 X 10 -^ 

7.0 X 10-2 

8.1 X 10-3 


0.005 

2.5 X 10 -^ 

8.3 X 10-2 

7.5 X 10-3 


0.002 

2.1 X 10 -^ 

1.1 X 10-3 

6.3 X 10-3 


0.1 

1.8 X 10-3 

5.3 X 10-2 

8.9 X 10-3 


0.05 

1.1 X 10-3 

6.1 X 10-2 

8.5 X 10-3 

LiteBIRD 

0.02 

7.3 X 10 -^ 

7.5 X 10-2 

8.0 X 10-3 

0.01 

5.8 X 10 -^ 

8.8 X 10-2 

7.1 X 10-3 


0.005 

4.9 X 10 -^ 

1.1 X 10-3 

6.2 X 10-3 


0.002 

4.4 X 10 -^ 

1.7 X 10-3 

5.3 X 10-3 


Table 4: la errors on r and rit for future satellite experiments. The pivot scale kp[Mpc 
is also listed correspondingly. 
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We plot the (r, rit) confidence ellipses of la and 2a C.L. for the CLASS and LiteBIRD 
experiments in Fig. Here the fiducial parameters are chosen as r = 0.01 and rit = 0. 
AdvACT and CLASS denote the best sensitivity on both r and rit in the five future 
ground-based experiments. Unfortunately, they can not still compare with the future 
satellite experiments such as LiteBIRD. Even though the satellite experiments give the 
best sensitivity on rit, it is impossible for them to test the consistent relation, i.e. rit = 
—r/8. By considering the BKP constraint ro,o 5 < 0.12 at 95% C.L., we deduce rit < 0.015 
which is much smaller than the la errors on rit in Tab. 



- 1 . 0 |- - 
0.000 0.005 0.010 0.015 0.020 

r 

Figure 1: The (r, rit) confidence ellipses of la and 2a C.L. for CLASS (blue) and LiteBIRD 
(red). The fiducial parameters are given by r = 0.01 and rit = 0. 


5 Discussion 

In this paper, we forecasted the sensitivity on the tensor tilt rit for future ground-based, 
balloon-borne and satellite experiments of B-mode polarizations. We focused on the scale- 
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invariant power spectrum = 0 which approximately corresponds to the canonical single- 
held slow-roll inhation models. For the tensor tilt, the satellite experiments can give a 
good sensitivity on rit even for r ~ 2 x 10“^. Specihcally, the la error on n* varies from 
~ 0.1 to ~ 0.01 when r varies from ~ 0.001 to ~ 0.1. By contrast, the ground-based and 
balloon-borne experiments give much worse sensitivity on both r and rit. Furthermore, 
our results did not change signihcantly with respect to previous forecasts on the future 
sensitivity of tensor-to-scalar ratio r in Ref. [25]. The reason is that we just considered the 
scale-invariant tensor tilt rij = 0 in our hducial model. By considering the BKP constraint 
on r, we see that it is impossible for these future experiments to test the consistent relation 
fit = —r/8 in the canonical single-held slow-roll inhation models. 
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A Instrumental specifications 


Experiments 

fskvl%] 

1/ [GHz] 

^FWHW ['] 

5P[^iK'] 

AdvACT 

50 

50 

90 

150 

2.2 

1.3 

7.8 

6.9 


50 

230 

0.9 

25 

CLASS 

70 

70 

38 

93 

90 

40 

39 

13 


70 

148 

24 

15 


70 

217 

18 

43 

Keck/BICEP3 

1 

1 

95 

150 

30 

30 

9.0 

2.3 


1 

220 

30 

10 

Simons Array 

20 

90 

5.2 

15.2 

20 

150 

3.5 

12.3 


20 

220 

2.7 

23.6 

SPT-3G 

6 

95 

1 

6.0 

6 

150 

1 

3.5 


6 

220 

1 

6.0 






EBEX 

1 

150 

8 

5.8 

1 

250 

8 

17 


1 

410 

8 

150 

Spider 

7.5 

94 

49 

17.8 

7.5 

150 

30 

13.6 


7.5 

280 

17 

52.6 







70 

30 

26 

19.2 

CMBPol 

70 

70 

45 

70 

17 

11 

8.3 

4.2 


70 

100 

8 

3.2 


70 

150 

5 

3.1 


70 

220 

3.5 

4.8 


70 

340 

2.3 

21.6 


70 

45 

23 

9.1 

COrE 

70 

75 

14 

4.7 

70 

105 

10 

4.6 


70 

135 

7.8 

4.6 


70 

165 

6.4 

4.6 


70 

195 

5.4 

4.5 


70 

225 

4.7 

4.6 


70 

255 

4.1 

10.5 


70 

285 

3.7 

17.4 


70 

315 

3.3 

46.6 


70 

375 

2.8 

119 


70 

60 

32 

10.3 

LiteBIRD 

70 

78 

58 

6.5 

70 

100 

45 

4.7 


70 

140 

32 

3.7 


70 

195 

24 

3.1 


70 

280 

16 

3.8 


Table 5: Instrumental specifications of future ground-based, balloon-borne and satellite 
experiments [32l [33l |3ll |35l |36l [371 EH EH HQl SI] • Here 5P = CTpixOpwHM- 
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